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SECTION  I 


INTRODUCTION 

Project  GEST  (Gas  Explosive  Simulation  Technique)  consisted  of  a  series 
of  tests  intended  to  simulate  late-time,  low  altitude  nuclear  fireball  phenome¬ 
nology.  Individual  events  involved  the  explosion  of  a  large  balloon  filled 
with  a  detonable  (near  stoichiometric)  mixture  of  methane  and  oxygen.  A  gaseous 
mixture  was  chosen  as  the  energy  source  rather  than  a  conventional  high  explo¬ 
sive  (TNT,  AN/FO,  etc)  for  two  reasons.  Gas  explosives  are  of  low  density  so 
that  there  is  no  problem  with  jetting  of  hot  species  out  of  the  fireball.  In 
addition,  gases  tend  to  form  a  more  stable  initial  fireball. 

In  a  nuclear  detonation  the  early  escape  of  energy  is  radiative.  As  the 
temperature  drops  below  a  few  thousand  degrees,  however,  the  energy  loss  by 
radiation  decreases  drastically.  The  hot  fireball  now  expends  energy  by  doing 
PdV  work  on  the  surrounding  atmosphere  and  grows  until  it  reaches  pressure 
equilibrium.  The  result  is  that  there  is  now  a  hot,  underdense,  mixture  of 
gases  and  bomb  materials  in  the  atmosphere.  Due  to  bouyancy  the  fireball  starts 
rising.  The  flow  around  the  rising  fireball  deforms  it  until  it  becomes  finally 
a  fully  developed  torus. 

In  rising  through  the  atmosphere  the  fireball  entrains  air  into  itself. 

This  entrained  air  interacts  with  the  bomb  debris  to  form  hot  metal  oxides 
which  radiate  in  the  10-12  p  wavelength  region.  Since  the  atmosphere  is 
transparent  to  these  wavelengths  the  effects  of  the  radiated  energy  are  felt 
at  a  considerable  distance  from  the  detonation. 

The  mechanisms  by  which  air  entrainment  and  subsequent  mixing  with  bomb 
debris  take  place  are  not  completely  understood.  Turbulence  at  the  edge  of  the 
fireball  and  subsequent  torus  may  be  a  primary  factor  in  the  entrainment  process. 
Knowledge  of  the  details  of  the  air  entrainment  and  mixing  processes  would  be 
of  great  value  in  determining  the  rate  of  metal  oxide  formation  and  the 
consequent  long  wavelength  infrared  (LWIR)  emission. 

Fireball  radar  clutter  is  another  facet  of  late  time  nuclear  fireball 
phenomenology  of  great  interest.  Its  source  mechanism  could  be  high 
temperatures  and  high  temperature  gradients  at  the  edge  of  the  torus,  or  it 
could  be  due  to  reflections  from  turbules  of  size  comparable  to  the  radar 
wavelength.  Although  it  was  not  expected  that  the  questions  of  radar  clutter 
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origin  could  be  fully  answered  in  the  types  of  experiment  discussed  here,  some 
additional  insight  could  perhaps  be  gained. 

The  role  of  turbulence  as  a  mechanism  in  the  various  phenomena  of  interest 
has  not  been  resolved.  Most  of  the  production-type  hydro  codes  are  laminar  flow 
codes  and  do  not  take  turbulence  into  account.  A  study  of  the  6EST  series  test 
results  in  comparison  with  HULL  predictions  should,  hopefully,  enable  a 
determination  as  to  the  desirability,  or  necessity,  for  introducing  turbulence 
into  the  codes. 

A  brief  description  of  the  experimental  setup  for  the  five  GEST  events 
involving  10  m  diameter  balloons  filled  with  a  detonable  mixture  of  methane 
and  oxygen  is  given  in  Section  II. 

Section  III  discusses  the  theoretical  modeling  of  the  GEST  event  using 
the  HULL  code.  Preliminary  results  from  the  diagnostic  instrumentation  which 
was  fielded  in  an  attempt  to  better  understand  the  fireball  and  torus  phenom¬ 
enology  are  included  in  Section  IV. 

The  last  section,  Section  V,  presents  the  preliminary  conclusions  which 
have  been  drawn  from  the  first  five  GEST  events.  Some  comparisons  are  made 
between  predicted  and  measured  values  of  torus  and  fireball  parameters.  These 
comparisons  reflect  favorably  on  the  HULL  code. 

A  sixth  GEST  event  -  a  double  burst  -  is  tentatively  scheduled  for 
February  1975.  This  will  be  a  repeat  of  the  fifth  GEST  event  which  was  a 
double  balloon  shot.  Because  of  operational  difficulties  one  of  the  two 
balloons  was  not  charged  with  the  correct  amount  of  methane  and  the  primary 
objective  of  that  particular  experiment  -  to  observe  air  flow  interactions  for 
multi  bursts  -  was  not  met. 

In  general,  the  first  four  GEST  events  which  involved  the  detonations  of 
single  balloons  can  be  considered  to  have  been  successfully  conducted  and 
much  useful  data  was  obtained. 
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SECTION  II 


THE  GEST  EXPERIMENT 

This  section  gives  only  a  brief  description  of  the  GEST  experiments  and 
related  technologies.  The  final  report  on  the  test  series  will  contain  more 
detailed  descriptions  of  the  physical  equipments  and  their  use  and  deployment. 

1.  TEST  SITE  LOCATION 

The  GEST  series  was  carried  out  at  the  Sandia  Laboratories  Aerial  Cable 
Test  Facility  in  the  southeast  field  test  areas  of  Kirtland  Air  Force  Base. 
Ground  zero  location  was  34°58'  North  latitude  and  106°27‘  West  longitude. 
Figure  1  shows  the  location  of  the  test  site  with  respect  to  topographical 
and  cultural  features  of  the  surrounding  area. 

2.  TEST  SITE  LAYOUT 

The  Sandia  Cable  Test  Facility  is  characterized  by  having  a  1-1/4  inch 
diameter  steel  cable  strung  across  a  relatively  flat  but  sloping  area  between 
hills  on  either  side  of  a  canyon.  The  maximum  height  to  which  the  cable  can  be 
raised  above  surface  ground  zero  (SGZ)  is  193  feet.  Figure  2  is  a  plan  view  of 
the  area  showing  the  locations  of  camera  stations,  instrumentation  vans,  etc. 
with  respect  to  the  cable  and  SGZ. 

3.  PHYSICAL  SETUP 
a.  Balloons 

The  envelopes  of  the  balloons  used  for  the  mair  GEST  events  were 
constructed  of  three  mil  Mylar.  The  Mylar  was  coated  on  both  sides  with  a 

O 

5000  A  thick  layer  of  aluminum.  They  were,  when  inflated,  10  meters  in 
diameter  and  had  a  volumetric  capacity  of  17,157  cubic  feet.  Figure  3  shows 
one  of  the  balloons  undergoing  leakage  testing. 

The  balloons  were  filled  with  a  mixture  of  methane  (CHi*)  and  oxygen 
(02)  in  the  approximate  ratio  of  1  part  to  2  parts,  respectively,  by  volume. 

The  detonation  of  this  mixture  in  the  quantity  assumed  was  approximately 
equivalent  in  energy  yield  to  the  detonation  of  1500  pounds  of  TNT. 

The  measured  quantity  of  oxygen  was  first  pumped  into  the  balloon 
through  a  fitting  and  poppet  valve  located  at  the  lowermost  point  of  the 
balloon.  After  the  required  amount  of  02  was  in  the  balloon,  the  methane  was 
pumped  in.  A  diffuser  installed  on  the  poppet  valve  acted  to  introduce 
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Figure  1.  GEST  Test  Site  Location 
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turbulent  motion  in  the  gas  flow.  The  lighter  CH4  bubbled  up  through  the 
heavier  02  in  a  light  plume  which  effectively  caused  mixing  of  the  gases  in  a 
reasonably  short  time  (reference  1). 

After  filling,  the  balloon  was  hoisted  to  the  detonation  altitude.  See 
table  1  which  lists  the  schedule  of  GEST  events,  the  relative  configurations, 
and  the  computed  yields.  Events  MB-1  through  MB-4  were  single  balloon  events; 
event  DMB-1  was  a  double  balloon  event. 

TABLE  1 .  SCHEDULE  OF  GEST  EVENTS 


EVENT 

DATE 

BURST  ALT. 
(FT) 

SEPARATION 

(FT) 

YIELD 

(109  JOULES) 

MB-1 

28  Nov 

73 

150 

NA 

2.86 

MB-2 

30  Nov 

73 

140 

NA 

2.93 

MB- 3 

11  Dec 

73 

130 

NA 

3.00 

MB-4 

13  Dec 

73 

130 

NA 

2.72 

DMB-1 

19  Jan 

74 

130 

100 

1.3-1. 6 
and  3.20 

Detonation  of  the  explosive  mixture  was  initiated  by  three  RP-1 
detonators  in  each  balloon.  Two  of  the  detonators  for  each  balloon  were  sur¬ 
rounded  by  plastic  sacks  each  containing  2.5  pounds  of  5-10  y  aluminum  powder. 
This  powder  was  dispersed  through  the  balloon  when  the  detonators  were  fired, 
and  acted  to  render  the  late  fireball  and  the  torus  optically  visible. 
Subsequent  to  the  second  shot  (MB-2)  50  gm  of  C-4  was  wrapped  around  each  of 
the  powder  dispersing  detonators  to  improve  powder  dispersal. 

4.  INSTRUMENTATION 

The  Air  Force  Weapons  Laboratory  (DEX)  fielded  time  of  arrival  (TOA) 
gages  and  blast  pressure  gages  on  each  of  the  five  events.  Experiments 
involving  IR  measurements  and  hot-wire  anemometry  were  fielded  by  AFCRL,  GE, 
LASL  and  TRW.  These  experiments  are  described  in  Section  IV,  together  with 
some  of  the  test  results  and  early  conclusions.  A  laser  experiment  fielded  by 
AFWL,  but  aborted  due  to  equipment  malfunction,  to  measure  Raman  scattering 
from  the  fireball  and  torus  is  also  described  there. 
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5.  PHOTO  COVERAGE 

Figure  2  gives  the  location  of  the  various  camera  stations.  Table  2  lists 
the  various  types  of  photo  coverage  at  each  of  the  stations.  Good  photo  footage 
was  obtained  on  all  events  and  a  documentary  film  is  in  preparation  as  part  of 
the  final  report.  Figures  4  thru  9  are  a  sequence  of  photos  of  the  MEM  event 
from  prior  to  detonation  to  formation  of  a  "picture  book"  torus. 
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Table  2.  Photographic  Coverage  of  6EST  Events 
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6b  MB-4,  DMB-1  70  ran  Multidata  10  5"  400'  B&W  T=3  min. 

(Dist=3808.8' )  DMB-1  16  ran  Milliken  32  15  inn  200*  color  T=5  min 
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Figure  4.  MB-4  Balloon  Prior  to  Detonation 
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Figure  (5.  MB-4  4  Seconds  after  Detonation 
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Figure  7.  MB-4  8  Seconds  after  Detonation 
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Figure  8.  MB4  12  Seconds  after  Detonation 
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Figure  9.  MB-4  Torus  60  Seconds  after  Detonation 
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SECTION  III 


CALCULATIONS 


1.  EARLY  CALCULATIONS 

Prior  to  the  GEST  series  a  computer  simulation  of  the  events  was  performed 
using  the  HULL  code.  These  calculations  are  described  in  reference  2.  Briefly, 
however,  the  calculations  involved  the  following.  The  SAP  code  (reference  3) 
was  first  used  to  burn  a  methane  oxygen  mixture  with  the  burn  being  "forced"  in 
that  the  detonation  wave  front  was  advanced  at  the  experimentally  determined 
wave  velocity  of  2.75  x  10s  cm/sec.  The  equation  of  state  used  for  the  unburned 
methane  oxygen  mixture  was  that  of  an  ideal  gas  with  y  =  1.3986  and  for  the 
burned  gas  y  was  taken  to  be  1.214.  For  the  equation  of  state  of  the  surround¬ 
ing  atmosphere  the  Doan-Nickel  equation  (reference  4)  was  used.  The  burn 
calculation  was  carried  out  until  the  resulting  shock  was  well  into  the 
surrounding  air. 

The  SAP  burn  output  was  then  used  as  initial  conditions  for  the  HULL  code 
portion  of  the  calculation.  The  SAP  data  were  inserted  into  an  ambient 
atmosphere  modeled  by  the  U.  S.  Standard  Atmosphere  with  the  balloon  center 
being  placed  70  m  above  the  ground  in  the  2-D  cylindrically  symmetric  HULL  mesh. 

The  chemistry  describing  the  detonation  was  expressed  using  the  MACE  code 
(reference  5).  With  thermodynamic  parameters  of  the  chemical  species  of  the 
problem  as  inputs  the  equilibrium  concentrations  were  calculated  by  specifying 
the  pressure  and  temperature  representative  of  the  specific  problem. 

2.  NEW  CALCULATIONS 

The  HULL  calculation  reported  here  is  similar  to  that  of  reference  3 
except  for  the  height  of  burst  and  the  equation  of  state.  For  the  new 
calculations  the  detonation  velocity  was  taken  to  be  2.41  x  10s  cm/sec  in  the 
SAP  burn.  The  actual  site  elevation  of  about  1800  m  above  sea  level  was  used  to 
determine  the  atmospheric  conditions  for  the  SAP  calculation.  These  conditions 
were  also  used  to  initialize  the  HULL  calculations.  This  new  HULL  calculation 
was  placed  43  meters  above  a  reflecting  surface  rather  than  70  m  as  was  the 
old  calculation.  Table  3  compares  the  old  and  new  sets  of  calculations. 
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TABLE  3.  COMPARISON  OF  OLD  AND  NEW  GEST  CALCULATIONS 


Original 

Calculations 

Final 

Calculations 

Height  of  Burst 

70  m 

43  m 

Atmospheric  Conditions 

STP 

Site  Conditions 

SAP  Calculation 

Burned  gas  equation  of  state 

Y=Const. 

y=y(p»e) 

HULL  Calculation 

Equation  of  State 

Air 

Air/Burned  methane 

In  order  to  incorporate  the  equation  of  state  for  burned  methane  in  the 
HULL  calculations,  it  was  necessary  to  identify  those  portions  of  the  fluid 
which  were  air  (material  1)  and  those  which  were  burned  methane  (material  2). 
This  multimaterial  calculation  was  run  with  a  "diffusion  limiter"  which  attempts 
to  preserve  material  gradients.  Although  these  procedures  can  give  rise  to 
spurious  flows  at  material  interfaces,  they  do  not  affect  the  total  flow  in 
the  calculation  since  only  a  few  cells  are  affected.  The  handling  of  two 
materials  in  this  fashion  affords  the  opportunity  to  determine  where  the 
detonation  products  are  as  a  function  of  time.  Figures  10,  11,  12,  and  13  are 
contours  in  the  r,z-plane  of  the  methane  burn  products  at  .5,  1,  2,  and  3 
seconds,  respectively. 

The  importance  of  the  lower  (43  m)  height  of  burst  is  illustrated  by 
figure  14.  This  is  a  plot  of  rise  velocity  during  the  first  second.  The 
large  rise  and  fall  in  velocity  at  .2  seconds  results  from  the  passage  of  the 
ground  reflected  shock  through  the  fireball.  Since  the  shock  is  stronger  for 
this  lower  (43  m)  HOB  than  it  was  for  the  older  calculation  (70  m  HOB),  the 
cumulative  rise  and  earlier  torusing  are  rather  significant. 

Figures  15,  16,  17  are  plots  of  the  altitude,  radius  and  velocity  of 
the  fireball  over  the  first  100  seconds.  Since  no  wind  profiles  are  included 
in  the  calculation,  a  large  deviation  in  rise  at  late  times  can  be  expected. 

This  matter  will  be  addressed  when  the  final  data  reductions  are  completed. 

The  most  dramatic  effect  of  using  the  methane  burn  products  equation  of 
state  is  reflected  in  the  temperature  time  history  (figure  18).  The  behavior 
of  the  burn  products  is  quite  different  from  that  of  air  at  a  given  internal 
energy  and  density.  The  new  calculation  shows  a  peak  temperature  of  about 
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Figure  10.  Burned  Methane  Density  Contours  at  0. 500  Seconds-HULL 
Calculation 
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RADIUS  M 

FIFWL  HULL  CALCULATION  OF  5  M  CH4  BALLOON  FIT  43  M  SAP  INPUT 
TIME  1.000  SEC  CYCLE  2142.  PROBLEM  1.0018 

Figure  11.  Burned  Methane  Density  Contours  at  1. 000  Seconds-HULL 
Calculation 
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Figure  12.  Burned  Methane  Density  Contours  at  2. 000  Seconds-HULL 
Calculation 
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Figure  13.  Burned  Methane  Density  Contours  at  3.000  Seconds-HULL 
Calculation 
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Figure  14a.  Early  Time  Rise  Velocity  and  Altitude  Versus  Time  -  43m 
HULL  Calculation. 
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ALTITUDE  (M) 


VELOCITY  (M/SEC) 


MASS  WEIGHTED  RISE  VELOCITY 


AFWL  HULL  CALCULATION  OF  5M  CH4  BALLOON  AT  70M  SAP  INPUT  19  OCT  72 

Figure  14b.  Early  Time  Rise  Velocity  and  Altitude  Versus  Time  -  70m 
HULL  Calculation. 
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Figure  15.  Late  Time  Altitude  versus  Time-HULL  Calculation 
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Figure  16.  Torus  Radius  versus  Time-HULL  Calculation 
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Figure  17.  Late  Time  RiBe  Velocity  versus  Time-HULL  Calculation 
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Figure  18.  Peak  Temperatures  versus  Time-Old  and  New  HULL  Calculations 
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3100°  at  burn  time  and  a  corresponding  slower  cooling  in  time.  An  important 
point  is  that  the  new  temperature  calculations  are  representative  of  the  true 
temperatures,  whereas  the  earlier  ones  were  not.  The  effect  of  the  methane 
equation  of  state  on  temperature  gradients  (figure  19)  is  not  as  marked  except 
that,  as  one  would  expect,  the  gradients  persist  at  somewhat  higher  magnitudes 
for  longer  periods  of  time. 

Additional  calculations  of  single  material  (air)  at  the  43  m  HOB  and  2 
material  (air  &  burned  methane)  problems  are  being  done  to  evaluate  the  shock 
torusing  effects  and  to  provide  more  detail  in  the  fireball  interior  for 
comparisons  with  the  data.  These  will  be  reported  separately. 

In  the  new  calculations  the  1-0  equilibrium  code  ODE  (reference  6)  was 
used  to  obtain  equilibrium  compositions.  This  code  uses  a  Newton -Raphson 
iteration  procedure  to  solve  the  equations  of  conservation  of  mass  and  energy 
subject  to  the  relationships  describing  the  state  of  chemical  equilibrium  and 
minimization  of  total  system  free  energy.  The  chemical  system  involved  in  the 
calculations  may  be  any  combination  of  ideal  gases,  »olids  or  liquids. 

The  ratio  of  specific  heats,  y»  for  the  equation  of  state,  and  the 
specific  heat  at  constant  volume,  Cy,  were  obtained  as  functions  of  density 
(p)  and  energy  (E).  To  facilitate  running  the  subsequent  SAP  and  HULL 
calculations  the  non-linear  least  squares  code  VARMIT  (reference  7)  was  used  to 
fit  a  surface  to  the  data.  The  surface  was  expressible  in  the  form 


Figures  20  &  21  show  constant  temperature  and  constant  temperature 
gradient  contours,  respectively,  for  T  +  2.0  seconds  as  calculated  with  the 
HULL  code  using  the  two  material  equation  of  state.  Superimposed  on  each  of 
these  figures  are  the  hot  wire  anemometer  strut  locations.  These  locations 
have  been  adjusted  vertically  to  reflect  the  fact  that  the  actual  height  of 
burst  for  event  MB-4  was  39  m  rather  than  the  43  m  for  which  the  computations 
were  run.  Figure  22  is  Included  here  to  illustrate  the  considerably  different 
temperature  contours  which  results  when  the  equation  of  state  for  air  (with 
constant  ratio  of  heat  capacities  =  1.4)  and  HOB  of  70  m  are  used  in  the  HULL 
computations. 
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TEMPERATURE  GRADIENT  (*  K/CM) 


Figure  19.  Temperature  Gradients  versus  Time-Old  and  New  HULL  Calculations 
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AFWL  HULL  CALCULATION  OF  5  M  CH4  BALLOON  AT  43  M  SAP  INPUT 
TIME  2. 000  SEC  CYCLE  4096.  PROBLEM  1. 0018 

Figure  20.  Constant  Temperature  Contours  at  2.000  Seconds-HULL 
Calculations 
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AFWL  HULL  CALCULATION  OF  5  M  CH4  BALLOON  AT  43  M  SAP  INPUT 
TIME  2.000  SEC  CYCLE  4096.  PROBLEM  1.0018 

Figure  21.  Constant  Temperature  Gradient  Contours  at  2. 000  Seconds-HULL 
Calculations 
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Figure  22.  Constant  Temperature  Contours  at  2. 000  Seconds-Old  HULL 
Calculations 
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SECTION  IV 


INDIVIDUAL  EXPERIMENTS 

In  addition  to  the  air  blast  pressure  and  time  of  arrival  measurements 
carried  out  by  AFWL/DEX  several  experimental  measurement  programs  were  carried 
out  by  groups  external  to  the  AFWL.  These  programs  used  or  intended  to  use 
infrared  (IR)  diagnostics,  and  hot-wire  anemometry  to  examine  fireball 
temperatures  and  species  concentrations,  turbulence  effects  on  species 
concentrations  and  on  temperatures  and  velocities.  These  individual  experiments 
and  certain  preliminary  findings  are  presented  below.  Where  they  were  available, 
preliminary  reports  by  the  individual  investigators  serve  as  basis  for  the 
material  presented. 
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1 .  TEMPERATURE  MEASUREMENTS 

(Los  Alamos  Scientific  Laboratories  -  Reference  13) 

IR  measurements  of  the  fireball  and  torus  were  made  by  LASL  during 
several  of  the  GEST  events.  Fireball  temperatures  could  be  deduced  from 
these  measurements.  Figures  23  and  24  give  the  temperatures  determined  for 
events  MB-1  and  MB-2.  The  device  used  was  an  ASA  thermovision  camera.  The 
camera  scans  the  field  of  view  and  focuses  the  infrared  radiation  on  a  detector 
which  converts  the  IR  signal  to  an  electrical  signal.  This  electrical  signal 
is  then  amplified  and  used  to  control  the  electron  beam  of  a  TV  monitor  tube, 
thereby  showing  the  termal  distribution. 

SAGA  is  a  short  wavelength  AGA.  It  uses  an  indium-antimonide  detector 
with  sensitivity  in  the  4.5  -  5.5  y  region.  LAGA  is  a  long  wavelength  AGA 
with  its  main  sensitivity  in  the  9.2  -  11.7  y  region;  however,  it  also  has 
windows  at  3.5  and  2.1  y.  LAGA  uses  a  mercury- cadmium- tell uride  detector. 

The  background  brightness  temperature  lines  were  derived  by  assuming  that 
the  blank  sky  which  served  as  backdrop  for  the  GEST  events  was  a  solid  black- 
body  radiator.  From  the  AGA  device  output  when  pointed  at  the  sky  the  equiv¬ 
alent  blackbody  temperature  was  calculated. 

In  the  calculations,  allowance  was  made  for  atmosniic-'ic  attenuation. 
However,  no  allowance  was  made  for(l)  less  than  unity  emissivitiss ,  (2)  transient 
cloud  transparencies,  or  (3)  masking  of  the  hottest  regions  by  optically  thick 
outer  regions  of  the  cloud.  All  of  these  would  end  u;  raise  the  calculated 
temperatures. 
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PEAK  BRIGHTNESS  TEMPERATURE  (K) 


Figure  23.  Peak  Brightness  Temperature  versus  Time-LASL  AGA  Measurements- 
Event  MB-1 
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PEAK  BRIGHTNESS  TEMPERATURE  (K) 


2.  MEASUREMENTS  OF  INFRARED  SPECTRA  IN  THE  RANGE  1.56  to  4.67  MICRONS 

(Utah  State  University  and  Cambridge  Research  Laboratories  -  Reference  14) 

Infrared  spectral  measurements  of  the  five  GEST  events  were  made  with  a 
rapid  scanning  Michelson  interferometer  specially  constructed,  by  Block 
Engineering,  to  AFCRL  and  USU  specifications.  The  interferometer  was  coaligned 
with  both  a  super  8  mm  camera  and  a  fixed  spectral  bandwidth  radiometer,  hand 
tracked  on  the  target  through  a  sighting  telescope.  The  8  mm  camera  confirmed 
tracking  of  the  instruments,  while  the  radiometer  gave  an  absolute  calibration 
check  for  the  interferometer  and  also  provided  absolute  emission  levels  in  a 
specified  bandwidth.  Thermal  pulse  radiometers  were  also  used  for  the  final 
GEST  explosion  to  give  early  time  information.  The  ground  range  of  the 
instrumentation  was  3425  feet  from  ground  zero  in  all  cases. 

The  data  from  the  interferometer  is  summarized  in  this  report.  The 
typical  results  which  are  presented  are  only  a  small  portion  of  the  extensive 
data  obtained.  These  preliminary  versions  of  the  absolutely  calibrated 
Infrared  spectra  obtained  with  4  cm"1  resolution  in  the  2140  to  6400  cm-1  (4.67- 
1.56  microns)  region  are  subject  to  what  will  probably  be  small  corrections 
arising  from  more  detailed  calibrations  of  the  neutral  density  filters. 
Comparisons  (not  given  here)  between  the  interferometer  and  radiometer 
intensities  confirm  the  absolute  levels  given  in  the  interferometer  data. 

Each  spectrum  presented  is  the  result  of  a  single  scan  of  the  Interfero¬ 
meter.  Because  the  interferometer  operates  at  a  rate  of  14  scans  per  second, 
each  scan  reduces  to  a  spectrum  which  is  the  spectral  irradiance  integrated 
over  a  time  interval  of  approximately  0.07  second.  Since  the  instrument  field 
of  view  was  4.5°  (determined  by  its  50%  response  points),  all  the  radiation 
from  the  fireball  and  debris  cloud  was  included  in  the  measurements. 

In  the  data  that  follow,  the  calibrated  spectral  irradiance  at  the 
Interferometer  entrance  aperature  was  determined  in  absolute  units  (watts -cm"2 
•wavenumber"1)  by  employing  field  and  laboratory  calibrations  with  black  body 
sources.  In  addition,  for  some  of  the  data,  temperature,  pressure,  and  dew 
point  measurements  made  at  the  test  site  were  used  as  inputs  to  atmospheric 
transmission  models  (reference  8)  to  determine  the  test  configuration 
atmospheric  transmission.  This  transmission  was  used  to  convert  from  spectral 
irradiance  values  to  absolute  spectral  radiant  intensities  (watts •steradlan"1 
•  wavenumber"1)  of  the  source.  All  of  the  data  has  background  radiation 
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subtracted  out,  and  the  contribution  of  sunlight  reflected  from  the  fireball 
and  debris  cloud  appears  to  be  negligible. 

In  figure  25,  the  spectral  irradiance  of  the  double  balloon  explosion 
DMB-1  is  shown  at  0.2  seconds  after  ignition.  The  sharp  dip  at  wavenumbers 
between  2200  and  2400  cm'1  is  due  to  atmospheric  absorption  by  C02.  The  large 
dips  between  3500  and  4000  cm"1  are  due  to  absorption  by  atmospheric  water. 

The  irradiance  at  that  time  is  at  least  three  orders  of  magnitude  above 
background.  The  very  noisy  wings  at  wave  numbers  below  2200  cm"1  and  above 
6000  cm"1  arise  because  the  interferometer  spectral  response  had  dropped  to  such 
a  low  value  in  those  regions  that  the  signal  was  primarily  instrument  noise. 
Currently  the  spectrum  is  interpreted  to  be  that  of  radiation  from  hot  H20 
bands  centered  near  1600,  3700  and  6250  cm"1  (6.3,  2.7  and  1.6  microns). 

In  figure  26,  the  spectral  irradiance  from  DMB-1  is  shown  at  1.59  seconds 
after  ignition.  The  temperature  of  the  radiating  water  has  dropped  significant¬ 
ly,  causing  the  radiation  to  concentrate  at  the  band  centers.  In  the  2200  to 
2400  cm"1  region,  hot  C02  emission  from  the  2300  cm-1  (4.3  micron)  band  is  now 
observable  as  it  is  transmitted  through  the  sides  of  the  atmospheric  C02 
absorption  band. 

In  figure  27,  the  spectrum  from  DMB-1  is  shown  3.27  seconds  after 
ignition.  The  water  bands  have  coded  further  and  the  C02  radiation  is  now 
very  prominent,  decaying  at  a  much  slower  rate  than  that  from  the  water.  The 
C02  band  irradiance  peaks  are  off  scale  and  appear  as  a  weaker  second  trace  at 
half  scale  in  the  2140  to  2800  cm"1  region. 

The  spectrum  of  DMB-1  is  shown  at  7.54  seconds  after  ignition  in  figure 
28.  A  small  amount  of  radiation  from  the  3700  cm"1  (2.7  micron)  H20  band  is 
still  observable,  and  the  C02  band  is  still  very  prominent. 

The  overall  time  variation  of  the  spectral  irradiance  at  seven  selected 
wavelengths  in  the  spectrum  is  shown  in  figure  29.  To  aid  in  the  interpretation 
of  this  figure,  curves  have  been  sketched  through  four  of  the  seven  sets  of 
points.  At  times  where  points  are  missing,  the  data  Is  not  available  because 
selective  optical  filters  for  other  measurement  requirements  were  being  used. 

In  figures  30  and  31,  the  spectral  irradiance  observed  from  MB-3,  the 
third  single  balloon  explosion,  is  shown  at  0.11  and  5.15  seconds  after 
ignition.  These  spectra  are  very  similar  to  those  already  described  from  DMB-1. 
In  fact,  all  of  the  methane-oxygen  balloon  explosions  exhibited  similar  behavior. 


44 


t .  - 


45 


Figure  25.  Sjpectrai  Irradiance  for  Event  DMB-1  at  0. 20  Seconds 


Figure  26.  Spectral  Irradiance  for  Event  DMB-1  at  1. 59  Seconds 
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Figure  27.  ^ectral  Irradiance  for  Event  DMB-1  at  3.27  Seconds 
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The  time  variations  of  the  same  seven  spectral  elements  displayed  in 
figure  29  for  DMB-1  are  shown  in  figure  32  for  MB-3.  Again,  the  behaviors  of 
DMB-1  and  MB-3  irradiances  appear  very  similar.  Although  DMB-1  was  a  double 
shot,  its  irradiance  was  not  significantly  greater  than  that  observed  from 
MB-3  because  of  atmospheric  attenuation  in  the  water  bands.  This  was  due  to 
the  fact  that  there  was  twice  as  much  water  in  the  atmosphere  when  DMB-1  was 
fired  (water  density  was  7.3  g/m  at  830  mb  and  4.0°C  for  DMB-1  as  compared  to 
3.4  g/m  at  830  mb  and  5.0°C  for  MB-3).  Also,  since  one  of  the  balloons  on 
DMB-1  was  not  completely  filled  with  CH4,  the  explosion  produced  a  smaller 
signature  than  had  been  anticipated. 

For  all  events,  the  C02  radiation  (solid  line  in  figures  29  and  32)  was 
visible  for  the  longest  time.  Water  emissions  at  wave  numbers  nearest  the 
band  heads  (short  dashed  lines)  persisted  significantly  longer  than  those  at 
wave  numbers  far  from  the  band  heads  (long  dashed  lines).  This  is  to  be 
expected  since  the  source  temperature  falls  rapidly  as  a  function  of  time. 

The  atmospheric  transmission  models  of  McClatchey  and  Selby  (reference  8) 
were  used  to  generate  the  spectral  transmission  data  shown  in  figure  33  for 
the  path  between  the  explosion  and  the  interferometer  on  MB-3.  The  data  was 
then  convolved  with  the  interferometer  instrument  line  shape  to  give  a  spectrum 
directly  comparable  to  the  spectra  measured  by  the  interferometer.  The  trans¬ 
mission  spectrum,  along  with  the  range  data,  was  used  to  determine  the  spectral 
radiant  intensity  of  event  MB-3  at  0.11  and  5.15  second-  and  the  results  are 
shown  in  figures  34  and  35,  respectively.  At  the  absorption  band  centers, 
where  the  atmospheric  transmission  is  very  small,  the  radiant  intensity  has 
great  uncertainty.  For  this  reason  the  2300  cm"1  (4.3  micron)  and  3700  cm'1 
(.27  micron)  regions  have  been  masked  out. 

The  relative  behaviors  of  the  spectral  irradiances  of  the  five  GEST  events 
at  3414  cm"1,  a  strong  H20  emission  region,  are  shown  in  figure  36.  All  of  the 
irradiances  at  this  wave  number  are  similar  up  to  one  second.  (It  should  be 
remembered  that  corrections  for  atmospheric  transmission  would  increase  the 
relative  levels  for  DMB-1  to  a  little  above  those  for  the  single  balloon 
explosions.)  Considerable  variations  occur  after  1  second.  For  instances, 
the  MB-4  irradiance  decayed  much  more  slowly  than  that  from  the  other  event; 
probably  because  this  shot  occurred  in  very  calm  conditions.  It  is  possible 
that  the  variations  and  different  decay  rates  were  due  to  the  effects  of  wind 
and  turbulence. 
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MB-3  ATMOSPHERIC  TRANSMISSION  (Model  at  4  CM  Resolution) 
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Figure  33.  Atmospheric  Spectral  Transmission  for  Event  MB-3 
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The  behavior  of  the  sharp  C02  spike  at  2396  cm-1  is  shown  in  figure  37 
for  all  events  except  MB-4,  for  which  it  was  not  measured.  Again,  all  of  the 
explosions  produced  similar  results  with  the  small  fluctuations  matching 
closely  those  shown  in  figure  36. 

It  is  felt  that  the  fireballs  and  debris  clouds  are  optically  thick 
during  most  of  the  observation  period.  Therefore,  the  data  correspond  to  the 
radiation  characteristics  of  the  outer  shell. 

The  CH4-O2  spectra  presented  here  appear  to  be  in  general  accord  with 
that  predicted  by  F.  N.  Alyea  (reference  9)  although  detailed  comparisons  have 
not  yet  been  made.  Water  and  carbon  dioxide  are  apparently  the  predjminant 
radiators. 

The  spectra  have  good  resolution,  which  will  allow  good  temperature 
information  to  be  determined.  This  will  provide  an  excellent  base  against 
which  model  predictions  of  fireball  radiation  can  be  compared.  The  spectral 
data  indicate  that  the  bulk  of  the  radiation  is  from  molecular  radiators 
(C02,  H20)  rather  than  from  a  blackbody-like  continuum.  Modeling  of  the 
spectral  radiation  as  originating  from  these  molecules  will  allow  the 
temperatures  of  the  C02  and  H20  to  be  determined. 


Figure  37.  Relative  Spectral  Irradiance  at  2396  cm-1  -  GEST  Events 


3.  RAMAN  SCATTERING  EXPERIMENT  TO  DETERMINE  AIR  ENTRAINMENT  RATES 

(Air  Force  Weapons  Laboratory  -  Reference  15) 

One  of  the  experimental  objectives  of  GEST  was  to  determine  the  rate 
of  air  entrainment  into  the  turbulent  mixture  of  hot  gases  formed  after  the 
explosion.  Laboratory  experiments  indicated  that  it  would  be  feasible  to 
determine  this  rate  by  measuring  Raman  scattering.  An  experiment  was  prepared 
for  the  double  balloon  event  DMB-1;  however,  difficulties  were  experienced  with 
the  equipment  prior  to  shot  time  and  the  experiment  was  aborted. 

Raman  scattering  is  a  form  of  inelastic  light  scattering.  When  a  photon 
of  frequency  v  interacts  with  a  molecule,  a  shift  in  frequency  (Av)  character¬ 
istic  of  the  scattering  molecule  can  result.  The  probability  of  the  interaction 
taking  place,  while  quite  small,  is  directly  proportional  to  the  number  of 
molecules  or  scattering  centers  available.  Thus,  for  a  given  initial  light 
intensity  of  known  wavelength  the  intensity  of  Raman  scattered  light  is  directly 
proportional  to  the  concentration  of  the  scattering  molecules. 

Based  on  this  effect  the  concentration  of  a  given  gas  or  gases  in  the 
turbulent  atmosphere  resulting  from  an  explosion  can  be  determined.  If  the 
concentration  is  determined  as  a  function  of  time  the  rate  of  air  entrainment 
can  be  deduced. 

In  order  to  demonstrate  feasibility,  a  high-power,  pulsed  ruby  laser 
(Space  Rays  Model  0A- 100)  was  used  in  a  laboratory  experiment.  Although  this 
particular  laser  could  not  be  pulsed  rapidly  enough  to  determine  air  entrainment 
rates  it  was  felt  that  a  satisfactory  analysis  of  a  single  pulse  would 
demonstrate  the  validity  of  the  approach.  The  experimental  plan  was: 

(1)  to  optimize  laser  operation  so  that  the  laser  consistently 
produced  a  single  pulse  of  the  order  of  109  watts. 

(2)  to  demonstrate,  in  the  laboratory,  an  ability  to  detect  the 
Raman  scattered  light. 

(3)  to  field  the  experiment  on  a  GEST  shot. 

The  optical  pickup  system  consisted  of  a  fl.5  cassegrainian  telescope, 
an  RCA  4526  (111  spectral  response)  photomultiplier  tube  (PMT),  a  Fluke  413B 
power  supply,  a  Tektronix  454  oscilloscope  and  optical  filters. 

The  optical  filters  were  designed  to  transmit  Raman  scattered  light  from 
02  and  C02  as  follows: 
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(i)  The  shifted  line  for  6943A  incident  light  is  at  7783A  for  02. 

The  filter  used  was  a  Heliotek  #8874-19  spiked  filter,  centered  at 

o  o 

7780A  with  a  79.5%  peak  transmission  and  a  68A  half  band  width. 

*  o  o 

(ii)  The  shifted  line  for  6943A  incident  light  is  at  7683A  for  C02. 

The  filter  used  was  a  Heliotek  #£.874-19  spiked  filter,  centered  at 

o  0 

7694A  with  a  73.5%  transmission  and  a  32A  half  band  width. 

(iii)  In  order  to  insure  suppression  of  wavelengths  outside  the  desired 
band  an  additional  0.8y  filter  was  used  in  combination  with  each 
of  the  above.  The  combined  filter  system  had  a  transmittance  of  20% 

o  o 

for  7683A  and  36%  for  7783A.  The  attenuation  outside  of  the  desired 
band  was  of  the  order  of  10"7  to  10"8. 

The  oscilloscope  monitoring  the  PMT  output  was  delay  triggered  using  the 
output  pulse  of  the  Pockel  cell  supply  of  the  0A-100  laser.  A  Pacific 
Photometric  Instrument  model  2A32  amplifier  with  a  Lambda  Model  LT2095  AM 
power  supply  was  added  to  the  system  when  it  was  discovered  that  system  rise 
time  was  too  slow.  The  amplifier  has  a  1  ns  rise  time  and  a  factor  of  10 
amplification.  With  this  modification,  the  form  of  the  detected  pulse  was 
similar  to  that  of  the  laser  pulse. 

With  the  detector  placed  at  various  distances  from  the  laser  beam  the 
Raman  signal  from  02  in  air  was  found  to  be  a  constant  seven  millivolts  for  a 
400  Mw  incident  pulse.  This  observation  was  consistent  with  the  low  divergence 
of  the  Cassegrainian  telescope  used  to  collect  the  signal.  Optical  filters  were 
interchanged  to  verify  that  it  was,  indeed,  the  Raman  and  not  the  Raleigh 
component  which  was  being  observed. 

The  system  was  fielded  on  DMB-1.  The  laser  was  placed  approximately  1000 
feet  from  one  of  the  balloons.  The  optical  pickup  system  was  placed  328  feet 
from  the  balloon  at  right  angles  to  the  laser  beam.  The  time  delay  adjustment 
was  accomplished  by  looking  at  the  integrated  signal  when  the  laser  beam  passed 
within  a  few  feet  of  the  telescope.  Then  with  the  amplifier  in  the  PMT  circuit, 
the  Raman  pulse  was  observed.  The  system  was  aligned  so  that  the  laser  beam 
would  intercept  the  hot  turbulent  gases  from  the  balloon. 

During  alignment,  a  malfunction  in  the  laser  power  generator  apparently 
caused  damage  to  the  Pockel  cell  driver  associated  with  the  0A-100.  Although  the 
generator  was  repaired  the  laser  would  not  then  function  correctly.  Attempts 
to  correct  the  difficulty  failed  and  the  experiment  was  aborted. 
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Although  external  difficulties  prevented  demonstration  on  the  GEST 
experiment,  laboratory  experiments  indicate  that  Raman  scattering  is  a 
potentially  successful  remote  analysis  tool.  Whether  Raman  scattering  can  be 
used  to  investigate  the  rate  of  air  entrainment  in  an  explosion  thus,  has  not 
yet  been  determined.  It  is  hoped  that  the  opportunity  to  perform  the  experiment 
with  another  GEST  type  event  will  arise,  and  the  results  will  be  favorable,  as 
expected. 


4.  FAST  RESPONSE  TEMPERATURE  MEASUREMENTS 

(TRW  Systems  Group  -  Reference  16) 

The  purpose  of  this  task  was  to  conduct  flow  field  temperature  measure¬ 
ments  within  the  rising  fireballs  of  the  GEST  methane-oxygen  detonations. 
Measurements  were  obtained  on  the  single  balloon  events  MB-1  through  MB-4  and 
on  the  double  balloon  event  DMB-1.  These  measurements  were  made  using  an  array 
of  fast  response  hot  wire  probes. 

Twelve  sensors  were  used  for  each  event,  nine  of  them  being  operated  as 
heat  transfer  gages  (constant  temperature)  and  three  as  constant  current 
resistance  thermometers  (T^).  The  two  types  of  hotwire  probes  used  are 
shown  in  figure  38.  In  all  cases  the  platinum-iridium  wire  sensors  were  0.00025 
inches  in  diameter.  In  the  constant,  temperature  mode  the  wire  sensors  were 
maintained  at  a  temperature  of  approximately  900°C  by  the  anemometer  feedback 
bridge.  All  data  were  recorded  on  a  VR-3700B  Magnetic  Tape  Recorder  in  the  FM 
mode  (DC  to  10  kHz). 

In  the  time  span  over  which  shock  front  temperatures  were  measured  (0-0.5 
seconds  after  time  zero)  all  of  the  sensors  performed  satisfactorily  on  events 
MB-1,  MB-3,  and  MB-4.  Nine  returnr-d  good  data  on  event  MB-2,  and  seven  on  DMB-1. 
Fireball  temperature  data  was  not  obtained  on  MB-1,  and  three  sensors  recorded 
the  temperature  of  the  fireball  leading  edge  only  on  MB-2.  The  entire  fireball 
time  span  was  covered  for  events  MB-3,  MB-4,  and  DMB-1.  For  this  data  interval 
(1.8  to  4.1  seconds  for  MB-3  and  MB-4  and  1.6  to  4.6  seconds  for  DMB-1)  nine, 
ten,  and  seven  of  the  sensors  performed  satisfactorily  for  the  three  events 
respectively.  The  portion  of  the  fireball  time  span  covered  on  MB-2  was  0.7 
to  1.2  seconds. 

Ten  of  the  sensors  were  mounted  on  a  triangular  framework  called  the  main 
probe  strut.  The  other  two  were  mounted  on  an  auxiliary  strut.  Figure  39  is 
a  photograph  of  the  main  and  auxiliary  struts  attached  to  the  overhead  cable. 

With  the  vertical  leg  of  the  main  strut  triangular  framework  positioned  at 
approximately  40  feet  from  the  fireball  centerline,  the  two  instrumented  legs 
measure  data  both  in  the  azimuthal  and  radial  directions.  The  use  of  seven 
and  three  probes,  spaced  18  inches  apart  on  the  aximuthal  and  radial  legs, 
respectively,  provided  results  spanning  the  expected  dimensions  of  interest. 

All  wire  sensors  were  aligned  in  a  direction  orthogonal  to  the  vertical  plane 
through  the  overhead  cable.  The  length  of  the  probe  stem  mounts  (10  inches) 
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Figure  38.  Hot  Wire  Probe  Configurations  for  GEST  Temperature  Measurements 
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Figure  39.  Probe  Strut  Installation  -  Event  MB-4 


and  the  orientation  angle  for  the  radial  sensors  (15°)  were  selected  to 
minimize  possible  interference  effects  from  the  supports.  Test  configuration 
drawings  for  the  single  and  double  balloon  GEST  events  showing  positions  of  the 
temperature  probes  are  given  in  figures  40  and  41. 

Typical  sensor  time  histories  for  the  three  radial  sensors  installed  on 
MB-4  are  presented  in  figure  42  and  43.  Responses  of  the  heat  transfer  gages 
(channels  1  and  3)  and  the  temperature  gages  (channel  2)  to  shock  front  motion 
and  fireball  are  distinctly  different.  This  difference  is  a  direct  result  of 
the  heat  loss  equation  for  the  constant  temperature  hot  wire  which,  for  a 
rising  fireball,  is  given  approximately  by 

e2 3 4  'v  q  'v  (  cT  +  •sTu  )  (T  -  T) 

W 

Here  e  is  the  bridge  output  voltage,  q  is  the  sensor  heat  transfer  rate,  c  is 
a  constant,  u  is  the  local  flow  velocity,  T  is  the  wire  temperature  (constant 

W 

=800°C),  and  T  is  local  flow  temperature.  Although  the  bridge  power  signal 
represents  the  combined  effects  of  both  the  velocity  and  temperature  fields 
within  the  fireball,  its  dependence  on  velocity  is  relatively  small,  both 
because  of  the  square  root  dependence  on  flow  velocity  and  because  the  heat 
transfer  constant  c  is  rather  large.  Velocities  within  the  fireball  can  be 
determined  approximately  from  optical  data  as  well  as  from  statistical  analysis 
of  turbulence  data.  As  a  result,  the  heat  transfer  measurements  can  be  reduced 
to  provide  relatively  accurate  temperature  data,  with  a  5  kHz  frequency  response 
as  obtained  from  the  anemometer  feedback  bridge.  Such  data  complements  the 
results  obtained  from  the  constant  current  sensors  which,  because  of  sensor 
thermal  response  characteristics,  are  limited  to  a  frequency  response  of  1  kHz. 

Additional  fireball  characteristics  noted  from  the  figures  were: 

(1)  the  fireball  was  turbulent  and  non-homogeneous  in  nature  at  early 
times 

(2)  temperature  fluctuations  were  of  the  same  order  as  the  "mean 
temperature  potential" 

(3)  channel  2  temperatures  approached  1000°C 

(4)  the  mean  temperature  gradient  at  the  fireball  leading  edge  (dT/dx) 
was  approximately  equal  to  15°C/cm.  This  is  similar  to  the  predicted 
value  (10  to  15°C/cm).  Instantaneous  temperature  gradients  at  eddy 
interfaces  were  as  high  as  200  to  300°C/cm 
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Figure  40.  Thermal  Sensor  Locations  -  Events  MB-1  through  MB-4 
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FIREBALL  Figure  43.  Expanded  Thermal  Time  Traces  -  Event  MB -4 


(5)  ambient  air  eddies  (pockets)  were  near  the  fireball  leading  edge 

(6)  there  was  general  agreement  among  records  as  to  the  large  eddy 
structure 

(7)  a  substantial  ambient  air  heat  transfer  rate  (large  stem  velocity 
after  fireball  passage)  was  indicated. 


5.  INFRARED  VIDICON  OBSERVATIONS 

(General  Electric  Company  -  Reference  17) 

The  General  Electric  Company  Space  Sciences  Laboratory,  under  sponsorship 
of  the  Defense  Advanced  Research  Projects  Agency  (DARPA),  deployed  two  Infrared 
vidicon  camera  systems  to  observe  GEST  events  MB-3  and  MB-4.  These  vidicon 
systems  were  equipped  with  narrow  band  spectral  bandpass  filters  centered  at 
3.3  microns  to  provide  infrared  diagnostic  data.  Purpose  of  the  observations 
was  to  aid  in  evaluation  of  the  flow  characteristics  of  fireballs.  Of  particular 
interest  were  the  characteristics  of  torus  formation,  including  formation  rate 
and  torus  size  and  rate  of  rise.  Temperature  profiles  across  the  torus  were 
also  of  interest. 

The  predicted  properties  of  the  flow  field  associated  with  the  debris  of 
the  balloon  explosion  were  obtained  from  HULL  hydrocode  computations.  These 
computations  are  described  in  reference  3.  The  major  products  of  the  explosion, 
C02  and  H20,  exhibit  strong  infrared  signatures  at  wavelengths  between  2.5  and 
3.5  microns.  In  order  to  obtain  spectral  radiant  intensity  predictions,  the 
temperature,  pressure,  and  concentration  profiles  provided  by  the  calculation 
were  used  to  solve  the  radiative  transfer  equation.  These  predictions,  for  the 
center  of  the  bright  regions,  are  shown  in  figures  44  and  45  for  2  seconds  and  6 
seconds  after  detonation. 

Two  infrared  vidicon  cameras  were  employed  for  the  measurements.  The 
cameras  had  identical  spectral  bandpass  characteristics,  but  had  predicted 
dynamic  intensity  ranges  which  were  arranged  in  such  a  way  as  to  maximize  the 
total  range  of  source  Intensities  which  could  be  observed. 

System  1,  the  "Bel voir"  camera,  was  borrowed  from  the  Army  Night  Vision 
Laboratory  at  Fort  Belvoir,  Virginia.  It  is  a  General  Electric  27934R 
reticulated  Ge:Cu  vidicon.  Optics  consisted  of  a  10-inch  focal  length  f/2  lens. 
The  cryogenic  system  was  a  pressure  Dewar  with  controlled  LN2  feed  to  the 
sensor.  The  unvignetted  field  of  view  was  4.1°  x  4.1°. 

System  2,  the  "NRL"  camera,  was  supplied  by  the  Naval  Research  Laboratory. 
It  is  a  General  Electric  Z7876  Ge:Cu  vidicon  with  a  4-inch  focal  length  f/1 
lens.  The  cryogenic  system  required  hand  filling  the  LN2  supply  to  the  vidicon 
tube.  Unvignetted  field  of  view  was  5.0°  x  5.0°. 
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Figure  44.  Predicted  Event  Temperature  Profile  at  T  =  2.0  Seconds 
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Figure  45.  Predicted  Event  Temperature  Profile  at  T=6. 0  Seconds 
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Both  systems  used  standard  30  frame  per  second  television  rasters.  Data 
was  recorded  on  Ampex  VR7000  video  recorders.  The  calibration  source  in  each 
case  was  a  Barnes  Engineering  Company  Model  LL-140T  Blackbody  with  a  12  x  12- 
inch  radiating  area. 

Both  cameras  had  identical  spectral  bandpass  filters  with  the  transmission 
curve  given  in  figure  46.  This  spectral  bandpass  was  determined  as  a  trade-off 
between  optical  thickness  effects,  atmospheric  transmission,  required  dynamic 
range  and  the  spectral  response  of  the  vidicon.  The  vidicon  spectral  response, 
from  reference  10,  r  shown  in  figure  47. 

An  observation  wavelength  of  3.3  microns  was  chosen  for  both  cameras  in 
order  to  observe  hot  H20  bands  while  minimizing  the  complicating  effects  of 
radiation  from  other  debris  species.  At  this  wavelength,  the  predicted  radiant 
intensities  (figures  44  and  45)  are  about  3  watts/cm2  ster  micron  at  two  and 
six  seconds  after  detonation  respectively.  The  requirement  for  the  overall 
dynamic  range  of  the  instrumentation  for  the  measurements  of  interest  was,  thus, 
about  3000.  Normally,  the  maximum  achievable  dynamic  range  of  the  cameras  is 
about  60  to  100. 

The  predicted  noise  equivalent  spectral  radiances  for  the  Bel  voir  and  NRL 
cameras  were  determined  to  be  3.6  x  10"5  and  2.45  x  10“5  watts/cm2 -ster-mi cron, 
respectively,  under  the  measurement  conditions  of  the  GEST  tests  (reference 
3).  These  levels  denote  the  minimum  signal  intensities  that  can  be  resolved. 
Thus,  while  the  cameras  had  sufficient  sensitivity  to  observe  the  event  at 
late  times  (around  6  seconds),  they  did  not  have  adequate  dynamic  range  to 
encompass  the  full  range  of  event  radiance  variation  from  early  to  late  times. 

To  extend  the  useful  dynamic  range  of  the  cameras,  each  was  fitted  with 
a  removable  aperture  stop  consisting  of  a  metal  plate  with  a  small  hole  drilled 
in  it.  These  aperture  stops  were  held  in  front  of  the  camera  optics  by  hand 
so  that  they  could  be  quickly  removed  when  the  event  spectral  radiance  fell 
beneath  the  sensitivity  limit  of  the  stopped-down  camera,  as  evidenced  by  a 
fading  of  the  TV  image.  The  Bel  voir  camera  stop  consisted  of  a  plate  with  a 
0.138-inch  diameter  hole,  which,  when  placed  over  the  five-inch  diameter  optics 
allowed  a  net  transmission  of  0.076%,  forming  an  f/ 72  optical  system.  Similarly, 
the  four- inch  NRL  lens  was  stopped  with  a  0.66- inch  aperture,  giving  a  trans¬ 
mission  of  2.7%,  or  an  f/6  optic.  In  this  manner,  the  camera  dynamic  ranges 
were  extended  to  cover  the  entire  region  of  interest,  as  shown  in  figure  48. 
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Figure  46.  Vidicon  Spectral  Filter  Transmission  Calibration  Curve 
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Figure  48.  Vidicon  Measurement  Dynamic  Range 
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The  event  spectral  radiances  were  expected  to  be  within  the  dynamic  range  of 
at  least  one  of  the  cameras  over  the  entire  period  of  observation. 

The  camera  systems  were  calibrated  using  the  vidicon  radiometric  transfer 
function 

S/N  =  C  J  ANx(c)Rl(X)  Ta(X)(.T1o(X)  T1  f(A)  dA 

where  S/N  =  ratio  of  peak  signal  to  peak-to-peak  noise  amplitude 

ANx(c)  =  differential  spectral  radiance  of  calibration  source 

above  that  of  adjacent  background  in  watt-cmz-ster-micron_1. 
At  the  source  temperatures  typically  used  in  these 
calibrations,  AN^cy*sN^cj ;  that  is,  the  background 

radiance  is  negligible. 

R 1  (x)  =  vidicon  relative  spectral  responsivity  normalized  to  the 

absolute  value  of  vidicon  responsivity  at  Argf,  the 
center  of  the  effective  system  passband. 

T  (X)  =  atmospheric  spectral  transmission  over  the  calibration 

a  C 

path. 

and 

T1q(x)  and  T 1 X)  =  relative  optics  and  filter  transmission  normalized  to 

their  absolute  values  at  A 

The  radiometric  calibration  process  consists  of  measuring  the  calibration 
constant  C  in  terms  of  observed  values  of  S/N,  given  values  of  AN^cj,  and  the 
calculated  value  of  the  integral  in  the  preceding  equation.  The  calibration 
was  performed  at  the  instrument  settings  at  which  the  data  was  obtained. 

A  posttest  calibration  was  also  performed,  which  consisted  of  observing 
the  radiance  of  a  12  x  12-inch  blackbody  source  at  a  distance  of  approximately 
50  feet  from  the  vidicon  cameras.  Five  source  temperatures  we1  e  used.  The 
calibration  integrals  were  evaluated  using  atmospheric  transmission  over  a  50- 
foot  horizontal  path  length  at  atmospheric  pressure  and  a  temperature  of  25°F. 
Altitude  was  taken  to  be  6000  feet  and  relative  humidity  25%.  A  light  haze 
(slight  scattering  loss)  was  assumed.  From  reference  11,  an  atmospheric 
transmission  (T  (A)J  of  0.995  at  X  =  3.3  microns  was  assumed. 

a  C 

The  results  from  posttest  calibration  of  the  Bel  voir  camera  are  shown 
in  figure  49.  This  calibration  differed  from  predicted  performance  and 
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Figure  49.  Belvoir  Camera  Post-Mission  Calibration 
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pretest  laboratory  calibration  in  that  the  observed  threshold  at  S/N  =  1  was 
approximately  three  times  that  previously  noted.  Dynamic  range  was  similarly 
decreased.  This  was  probably  due  to  increased  noise  in  the  video  recordings 
taken  at  the  field  site. 

The  NRL  camera  vidicon  tube  was  quite  old  and  began  to  "go  soft"  after 
arrival  at  the  field  site.  A  replacement  was  not  available  in  time  for  the 
tests.  As  a  result,  noise  equivalent  source  radiance  was  an  order  of  magnitude 
higher  than  expected  and  dynamic  range  was  drastically  reduced.  A  noise- 
equivalent  source  radiance  of  1.62  x  10-4  watt/cm2-  ster  was  determined  from 
the  posttest  calibration  data. 

The  best  data  was  obtained  by  the  Bel  voir  camera  on  MB-4.  Because  the 
absence  of  a  boresite  lamp  on  the  cable  made  optical  alignment  impossible,  no 
data  was  obtained  by  this  camera  on  event  MB-3.  On  MB-4,  the  Bel  voir  camera 
observed  the  event  from  detonation  to  about  T  +  9.5  seconds.  The  system  was 
saturated  over  part  of  the  middle  region  as  predicted  by  the  experiment  plan 
(reference  12). 

The  NRL  camera  obtained  5  seconds  of  data  on  both  MB-3  and  MB-4.  This 
camera  was  successfully  aligned  for  MB-3  even  without  the  boresite  lamp. 

Despite  the  vidicon  degradation  discussed  previously,  the  NRL  data  is  of 
adequate  quality  to  allow  Identification  of  source  characteristics  at  the 
threshold  intensity,  which  for  this  camera  was  about  equivalent  to  a  465°K 
blackbody. 

A  sample  data  frame  is  given  in  figure  50.  The  data  is  displayed  and 
analyzed  on  a  video  tape  display  station  which  permits  a  single  video  frame 
to  be  shown  on  a  TV  monitor,  or  C-scope,  and  also  permits  a  single  video 
amplitude  of  a  single  TV  scan  line  to  be  displayed  on  an  oscilloscope,  or  A- 
scope.  The  TV  scan  line  which  has  been  selected  for  display  on  the  A-scope  is 
indicated  by  a  white  marker  pulse  on  the  C-scope.  The  position  of  the  marker 
pulse  on  the  scan  line  is  shown  as  a  pulse  on  the  lower  trace  of  the  A-scope. 
The  vidicon  tube  in  the  Bel  voir  camera  had  a  major  blemish  on  its  retina,  which 
shows  up  as  constant  spot  to  the  right  of  the  frame  center  in  the  C-scope 
photograph.  The  sample  data  frame  of  figure  50  shows  the  rising  cloud  at  a 
time  when  part  of  the  signal  is  above  the  saturation  limit  of  the  vidicon.  At 
saturation,  the  vidicon  exhibits  a  reversal  effect,  so  that  the  saturated 
portion  of  the  image  appears  black  on  the  C-scope.  The  black  areas  on  the  C- 
scope  correspond  to  the  minima  in  the  A-scope  video  signal  level. 


81 


82 


%  - 


Figure  50.  Sample  Vidicon  Data  Frame 


Typical  data  frames  from  the  exploding  balloon  event  MB-4  are  shown  in 
figure  51.  The  first  frame,  corresponding  to  time  of  T  +  0.2  seconds,  shows 
the  cloud  entering  the  bottom  of  the  camera  field-of-view.  Because  the  aperture 
stop  was  in  place  at  detonation,  the  observed  intensity  is  actually  attenuated 
by  a  factor  of  7.6  x  10"4.  The  next  two  frames,  at  time  T  +  1.1  and  T  +  1.6 
seconds,  show  the  rise  of  the  roughly  spherical  cloud.  During  this  period 
the  radiant  intensity  of  the  cloud  is  falling  rapidly,  and  shortly  after  the 
T  +  1.6  second  frame  intensity  falls  below  the  threshold  of  the  vidicon. 

At  this  point  in  time,  the  aperture  stop  was  removed  from  the  camera  lens. 
The  results  are  shown  in  part  2  of  figure  51.  The  frame  at  T  +  2.8  seconds  and 
the  next  frame  are  interesting  in  that  they  show  the  cloud  while  the  aperture 
stop  was  being  removed.  The  dynamic  range  distribution  shown  in  figure  48 
indicates  that,  for  the  Belvoir  camera,  there  is  an  interval  of  time  in  which 
the  radiant  intensity  is  below  threshold  for  the  stopped  vidicon  but  above  the 
saturation  limit  for  the  unstopped  vidicon.  The  aperture  stop  was  removed 
manually,  and  therefore  complete  removal  required  a  finite  time.  These  two  data 
frames  were  obtained  when  the  aperture  stop  was  only  partially  removed.  The 
cloud  radiant  intensity,  however,  falls  within  the  dynamic  range  of  the 
partially  stopped  vidicon.  These  frames  clearly  show  the  flattening  of  the 
cloud  and  the  beginning  of  torus  formation. 

When  the  aperture  stop  is  completely  removed,  at  T  +  3.1  seconds, 
saturation  obscures  the  details  of  the  central  cloud  structure. 

The  next  series  of  data  frames,  part  3  of  figure  51,  shows  further 
flattening  of  the  cloud  and  much  clearer  definition  of  the  torus.  By  T  +  5.0 
seconds,  the  torus  is  well  defined  with  a  very  conspicuous  size  and  intensity 
imbalance  between  its  right  and  left  sides.  The  early  stages  of  this  imbalance 
can  be  seen  in  the  T  +  2.8  frame. 

The  final  series  of  data  frames,  part  4  of  figure  51,  shows  the  fade-out 
of  the  fully  developed  torus.  The  frames  at  T  +  7.5  and  T  +  8.5  seconds  show 
the  fully  developed  torus  and  a  cross-section  of  the  much  brighter  left  side. 

A  feature  of  the  data,  not  apparent  in  the  still  frames  but  readily  discernable 
in  the  video  data,  is  a  rotating  or  "churning"  motion  about  the  circular  axis. 

At  T  +  9.5  seconds,  the  torus  simultaneously  rises  out  of  the  vidicon  field-of- 
view  and  falls  below  the  instrument  sensitivity  threshold. 
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Figure  51.  (Part  4)  MB-4  Vidicon  Data  Frames  (Continued) 
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Additional  analysis  of  the  data  from  MB-4  has  produced  the  results  shown 
in  figures  52  through  62.  Orientation  of  the  left  hand  toru:;  image  outline  at 
vidicon  threshold  as  a  function  of  time  is  given  in  figure  52.  Data  points 
from  both  the  Belvoir  (#1)  and  NRL  (#2)  cameras  are  included.  Data  from  the 
two  cameras  are  in  good  agreement  and  indicate  how  the  image  flattens  out  as 
the  torus  rises. 

figures  53  and  54  give  the  vertical  dimension  of  each  side  of  the  torus 
as  a  function  of  time.  The  image  outline  is  obtained  from  both  cameras  as  a 
contour  through  vidicon  threshold  intensity.  Imbalance  between  right  and  left 
sides  of  the  torus  is  very  apparent  in  figure  53.  Horizontal  torus  dimensions 
as  functions  of  time  are  given  in  figures  55  and  56.  Torus  separation,  total 
image  width,  and  cloud  mean  altitude  as  functions  of  time,  as  obtained  from 
the  Belvoir  camera  data,  are  given  in  figures  57,  58,  and  59.  Definitions  of 
all  of  these  quantities  may  be  obtained  from  the  small  diagrams  included  with 
each  figure. 

Spectral  radiant  intensities  across  the  horizontal  field  of  view  of  the 
camera  are  presented  in  figures  60  and  61.  Each  line  in  the  figures  represents 
a  single  TV  scan  line,  as  might  be  displayed  on  the  A-scope.  The  distance 
in  meters  of  each  scan  line  from  the  top  of  the  image  (threshold  intensity)  at 
the  image  plane  is  given  in  the  figure  key.  All  of  the  data  in  both  figures 
60  and  61  is  taken  from  a  frame  at  T  +  5.9  seconds. 

Finally,  the  iso-intensity  contour  at  2.5  x  10”3  watt/cm2-ster-micron  is 
plotted  in  figure  62.  This  contour  is  taken  directly  from  the  curves  of 
figures  60  and  61  at  T  +  5.9  seconds. 

The  data  obtained  by  the  Belvoir  camera  on  MB-4  is  of  excellent  quality 
and  in  general  the  experiment  is  considered  to  have  been  successful.  Additional 
analysis  is  planned  and  will  be  available  in  the  final  report.  A  16  mm  film 
of  the  vidicon  data  can  be  obtained  by  contacting  Rome  Air  Development  Center/ 
OCSE,  Attention:  Mr.  Joseph  J.  Simons,  Griffiss  Air  Force  Base,  Rome,  New  York 
13441  (Refer  to  Contract  F3602-73-C-0102) . 
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Figure  52.  Torus  Orientation  versus  Time  -  Event  MB-4 


Figure  53.  Torus  Region  Height  versus  Time  -  Event  MB-4 


Camera  1 


Figure  54.  Torus  Region  Height  versus  Time  -  Event  MB-4  -  Camera  2 
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Figure  55.  Torus  Region  Width  versus  Time  -  Event  MB-4  -  Camera  l 


92 
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Figure  56.  Torus  Region  Width  versus  Time  -  Event  MB-4  -  Camera  2 
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Figure  57.  Torus  Region  Separation  versus  Time  -  Event  MB-4  -  Camera  1 
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Figure  59.  Cloud  Mean  Altitude  versus  Time  -  Event  MB-4  -  Camera  1 
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AT  3.3  MICRONS  (WATT/CM* -STER-MICRON) 


Figure  60.  Spectral  Radiant  Intensity  versus  Horizontal  Distance  - 
Upper  Portion  of  Image  -  Event  MB-4 
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AT  3.3  MICRO!®  (WATT/CM2-STER-MICRON) 
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6.  AFWL  TEST  MEASUREMENT  RESULTS 

In  addition  to  the  contractor  conducted  experiments  described  above, 
measurements  of  air  blast  time  of  arrival  and  pressure  history  were  made  by 
AFWL/DEX. 

a.  Time  of  Arrival  Measurements 

For  each  of  the  MB  shots  22  PZT  crystals  were  fielded  to  obtain 
the  time  of  arrival  of  the  shock  front.  Four  of  these  crystals  were  placed 
in  pockets  on  the  balloon  envelope  and  six  were  mounted  on  each  of  three 
cables  as  shown  in  figure  63. 

Balloon  Time  of  Arrival  (BTOA)  gages  were  mounted  to  the  surface  of  the 
balloon  with  BTOA-1  being  3  feet,  BTOA-2  27  feet  6-3/4  inches,  and  BTOA-3  and  -4 
32  feet  0-3/8  inches  from  the  balloon  top  end  fitting.  The  gages  were  arranged 
to  be  at  the  apices  of  a  tetrahedron.  They  were  thus  5  meters  from  the  balloon 
center. 

Cable  Time  of  Arrival  (CTOA)  gages  were  cable  mounted  with  positions 
along  the  cables  so  determined  that  distances  from  the  gages  to  the  balloon 
center  were  as  follows: 


CT0A-X.1 

12  meters 

C\J 

X 

5 

14  meters 

-X.  3 

16  meters 

X 

1 

18  meters 

-X.5 

20  meters 

-X.6 

22  meters 

where  X  =  1,  2,  3  denotes  the  cable  number.  These  dimensions  are  valid  for  shot 
MB- 1  only  where  the  HOB  was  150  feet.  For  shot  MB-2  the  HOB  was  140  feet  and 
for  MB-3  and  -4  it  was  130  feet.  The  gages  were  mounted  at  the  same  locations 
on  the  cables  for  all  of  the  MB  series.  Figure  2  lists  the  locations  of  the 
ground  ends  of  the  TOA  cables  with  respect  to  S6Z. 

For  event  DMB-1  four  crystals  were  mounted  on  each  of  the  two  balloons. 
There  were  no  cable  mounted  transducers  for  this  event. 

Table  4  lists  the  measured  times  of  arrival  of  the  shock  front  at  each  of 
the  TOA  crystals.  The  results  of  these  measurements  were,  in  general,  not 
very  satisfactory  because  of  the  choice  of  sensor  (PZT  ceramic  crystal) 
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Figure  63.  GEST  TOA  Crystal  Deployment  Configuration 
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TABLE  4.  GEST  TOA  DATA 


EVENT 


XTAL 

MB-1 

MB- 2 

MB- 3 

M3-4 

DMB-1 

.1  DMB-1. 2 

BTOA-1 

1.993 

- 

1.773 

1.968 

- 

- 

-2 

1.938 

- 

- 

1.996 

- 

1.969 

-3 

1.912 

- 

- 

- 

2.175 

- 

-4 

- 

1.890 

2.098 

1.997 

- 

1.938 

CTOA-1.1 

- 

- 

- 

- 

NOT 

INSTRUMENTED 

-1.2 

- 

- 

16.951 

17.175 

-1.3 

- 

- 

19.035 

19.460 

-1.4 

- 

- 

21 .979 

- 

-1.5 

- 

- 

25.245 

25.358 

-1.6 

- 

- 

29.096 

29.071 

CTOA-2.1 

- 

- 

- 

- 

CVI 

CM 

1 

- 

16.130 

14.685 

13.631 

-2.3 

- 

19.130 

17.046 

15.761 

-.24 

- 

22.726 

19.691 

18.280 

-2.6 

- 

26.410 

22.868 

24.322 

-2.6 

- 

30.762 

26.743 

25.134 

CTOA-3.1 

9.55 

- 

- 

- 

-3.2 

13.78 

- 

14.691 

- 

-3.3 

- 

- 

17.525 

- 

-3.4 

- 

- 

20.883 

- 

-3.5 

- 

- 

24.466 

- 

-3.6 

_ 

_ 

_ 

NOT 

INSTRUMENTED 

102 


and  its  mounting  arrangement.  Crystal  outputs  tended  to  be  lower  than  desirable 
for  accurate  measurement.  The  times  presented  are  time  intervals  between  FIDU 
pulse  and  TOA  pulse  in  milliseconds. 

b.  Air  Blast  Pressure  Measurements 

Ten  Kulite  XT  Series,  integrated-sensor,  Q  type  pressure  transducers 
were  fielded  on  each  of  the  MB  events.  Nine  of  these  were  mounted  on  the  TOA 
cables  as  shown  in  figure  64;  one  was  located  at  surface  ground  zero. 

For  the  DMB-1  event  twelve  pressure  transducers  were  mounted  at 
ground  level  in  the  configuration  shown  in  figure  65.  Pressure- time  histories 
from  these  gages  are  shown  in  figures  66  through  77. 

At  this  writing  the  pressure  measurement  data  has  not  been  analyzed. 
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Figure  64.  GEST  Pressure  Gage  Deployment  Configuration  -  MB  Series 
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MEASUREMENT  NO. 

RANGE  (FT) 

GAGE  RANGE  (PSI) 

100 

0 

100 
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50 

50 

102 

50 

50 

103 

50  _L 

10 

104 

25J_ 

100 

105 

10J_ 

100 

106 

25 

50 

107 

10 

100 

108 

10 

100 

109 

25 

50 

110 

100 

10 

111 

200 

10 

Figure  65.  GEST  Evenl  DMB-1  Blast  Pressure  Gage  Layout 


105 


17.5 


(soas-isd)  asindrn 


©  m  ©  to  ©  »o 

•  •  •  •  •  *  •  •  • 

lO  Cl  ©  C**  C'j©c',u^ 

rH  iH  «H  11 


(isd)  annssand 


106 


Figure  66.  GEST  DMB-1  Pressure  Measurement  100 
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Figure  71.  GEST  DMB-1  Pressure  Measurement  105 
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Figure  72.  GEST  DMB-1  Pressure  Measurement  106 
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Figure  75.  GEST  DMB-1  Pressure  Measurement  109 
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Figure  76.  GEST  DMB-1  Pressure  Measurement  110 
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SECTION  V 


DISCUSSION  OF  RESULTS 

Generally  good  results  were  obtained  with  the  GEST  series  of  balloon 
explosions.  There  were  the  usual  experimental  difficulties  such  as  loss  of 
the  Raman  scattering  experiment  and  the  failure  to  place  the  required  quantity 
of  methane  in  one  of  the  DMB-1  balloons.  This  latter  failure  precluded 
verification  of  the  HULL  prediction  that  the  two  tori  would  interact  to  form 
a  single  larger  one. 

Hot  wire  anemometry  data  obtained  in  the  test  series  can  be  used  to 
obtain  temperature,  temperature  gradients,  temperature  fluctuations,  velocity, 
velocity  gradients,  and  velocity  fluctuations  for  comparison  with  HULL  calcu¬ 
lations.  From  the  fluctuation  data  turbulence  information  can  be  derived.  An 
approximately  value  for  the  ratio  Tf-j^/T  n  was  determined  for  the  GEST  shots 
to  be  0.3.  This  is  important  in  that  it  implies  that  the  temperature  gradients 
calculated  by  a  laminar  code  such  as  HULL  would  not  be  significantly  reduced  by 
the  inclusion  of  turbulence. 

Preliminary  study  of  the  turbulence  information  yields  an  estimate  of 
0.4  seconds  for  the  dissipation  time  and  about  4  meters  for  the  integral  scale. 
Although  preliminary,  these  figures  reflect  the  sort  of  information  obtainable 
from  GEST  type  events.  Figure  78  is  the  power  spectrum  derived  from  one  of 
the  TRW  hot  wire  measurements.  The  straightline  fit  over  much  of  the  frequency 
range  has  a  slope  of  -5/3.  That  this  slope  corresponds  to  that  expected  in 
the  inertial  subrange  for  isotropic  turbulence  is  reassuring  and  tends  to  rein¬ 
force  confidence  in  the  validity  of  the  data. 

The  TRW  data  showed  typical  values  of  15°  C/cm  for  the  temperature 
gradients  at  the  fireball  leading  edges  and  0.5°  C/cm  at  the  trailing  edges. 

The  GEST  type  detonation  with  its  high  temperature  gradients  appears  to  offer 
a  useful  tool  for  the  study  of  fireball  radar  clutter  and  its  origin. 

One  of  the  prime  reasons  for  carrying  out  the  GEST  series  of  shots  was 
that  it  would  give  data  which  could  be  used  to  verify  the  accuracy  of  the  HULL 
code  for  the  prediction  of  fireball  and  related  phenomenology.  Preliminary 
analysis  of  the  data  indicates  that  the  HULL  code  is  a  useful  tool  in  this 
regard.  Temperatures,  temperature  gradients,  torus  formation  rates,  etc.  as 
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predicted  using  the  43  m  HOB  and  the  2  material  equation  of  state  were  in 
reasonably  good  agreement  with  the  preliminary  data  determinations. 

Figure  79  compares  the  LASL  determined  temperatures  for  MB-2  using  the 
AGA  camera  with  the  HULL  predictions.  That  the  AGA  temperatures  are  low 
compared  to  prediction  is  probably  due  to  the  fact  that  LASL  made  no  allowance 
for  less  than  unity  emissivi ties,  transient  cloud  transparencies,  or  masking 
of  the  hottest  regions  by  optically  thick  outer  regions.  All  of  these  omissions 
are  in  the  direction  to  depress  the  measured  temperature  determination.  That 
both  the  5  p  and  10  p  measurements  are  so  close  to  each  other  is  probably  due  to 
the  sharp  temperature  gradients  being  observed. 

Figure  80  presents  a  comparison  of  HULL  predicted  torus  dimensions  as  a 
function  of  time  with  visual  (photographic)  data.  The  correspondence  is  very 
good. 

Figure  81  is  a  comparison  of  fireball  and  debris  cloud  altitude  versus 
time  as  observed  visually  (photographically)  and  as  calculated  by  the  HULL 
code.  Again,  the  correspondence  between  prediction  and  experimental  measure¬ 
ment  is  seen  to  be  quite  good. 

Torus  formation  on  all  single  shots  was  good  although  for  event  MB-2  the 
seedant  dispersal  was  quite  non-uniform  and  the  day  was  somewhat  cloudy  so  that 
satisfactory  optical  records  were  not  obtainable.  Infrared  observations, 
however,  confirmed  good  torus  formation. 


120 


EXTREME  AND  MASS  WEIGHTED  MEAN  RADIUS 


AFWL  HULL  CALCULATION  OF  5  M  CH4  BALLOON  AT  43  M  SAP  INPUT 


Figure  80  Torus  Radius  versus  Time-Optical  Measurements  and  HULL  Predictions 
Compared 
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Figure  81.  Altitude  versus  Time  -  MB  Series  and  HULL  Prediction  Compared 
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